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Abstract 

Estimating and controlling large risks has become one of the main 
concern of financial institutions. This requires the development of ad- 
equate statistical models and theoretical tools (which go beyond the 
traditionnal theories based on Gaussian statistics), and their practi- 
cal implementation. Here we describe three interrelated aspects of 
this program: we first give a brief survey of the peculiar statistical 
properties of the empirical price fluctuations. We then review how 
an option pricing theory consistent with these statistical features can 
be constructed, and compared with real market prices for options. 
We finally argue that a true 'microscopic' theory of price fluctuations 
(rather than a statistical model) would be most valuable for risk as- 
sessment. A simple Langevin-like equation is proposed, as a possible 
step in this direction. 

1 Introduction 

The efficiency of the theoretical tools devised to address the problems of 
risk control, portfolio selection and derivative pricing strongly depends on 
the adequacy of the stochastic model chosen to describe the market fluc- 
tuations. Historically, the idea that price changes could be modelled as a 
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Brownian motion dates back to Bachelier [T|]. This hypothesis, or some of 
its variants (such as the Geometrical Brownian motion, where the log of the 
price is a Brownian motion) is at the root of most of the modern results of 
mathematical finance, with Markowitz portfolio analysis, the Capital Asset 
Pricing Model (capm) and the Black-Scholes formula || standing out as 
paradigms. The reason for success is mainly due to the impressive mathe- 
matical and probabilistic apparatus available to deal with Brownian motion 
problems, in particular Ito's stochastic calculus. 

An important justification of the Brownian motion description lies in 
the Central Limit Theorem (clt), stating that the sum of N identically 
distributed, weakly dependent random changes is, for large N, a Gaussian 
variable. In physics or in finance, the number of elementary changes observed 
during a time interval t is given by N = 4r where r* is an elementary 
correlation time, below which changes of velocity (for the case of a Brownian 
particle) or changes of 'trend' (in the case of the stock prices) cannot occur. 
The use of the CLT to substantiate the use of Gaussian statistics in any 
caserequires that t ^> r*. In financial markets, r* turns out to be of the 
order of several minutes, which is not that small compared to the relevant 
time scales (days), in particular when one has to worry about the tails of the 
distribution (corresponding to large shocks) which sometimes disappear only 
very slowly. 

The fact that the tails of the distribution of returns are much 'fatter' than 
predicted by the gaussian is well known, in particular since the seminal work 
of Mandelbrot [|J, where the idea that price changes are still independent, 
but distributed according to a Levy stable law was first proposed. This model 
however fails in two respects. First, the tails of the distribution of returns 
is now much overestimated] in particular, the variance appears to be well 
defined in most financial markets, while it is infinite for Levy distributions. 
Second, and perhaps more importantly in view of its application to option 
markets, the amplitude of the fluctuations (measured, say, as the local vari- 
ance) appears to be itself a randomly fluctuating variable with rather long 
range correlations. 

The aim of this paper is to provide a short survey of the most promi- 
nent statistical properties of the fluctuations of rather liquid markets, which 
are characterized by what one could call 'moderate' fluctuations (for a re- 
view, see e.g. ||). Less liquid markets sometimes behave rather differently 
and more 'extreme' fluctuations can be observed. We shall then present a 
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theory for option pricing and hedging in the general case where the under- 
lying stock fluctuations are not gaussian. In this case, perfect hedging is in 
general impossible, but optimal strategies can be found (analytically or nu- 
merically) and the associated residual risk can be estimated. We show that 
the volatility 'smile' observed on option markets can be understood using a 
cumulant expansion, and discuss the idea of an implied 'kurtosis', which is 
(on liquid markets) very close to the actual (maturity dependent) kurtosis of 
the historical data. Finally, we briefly discuss a simple 'Langevin' approach 
to market fluctuations, which aims at describing, in a coarse-grained way, 
the feedback effects between the agents behaviour and the price fluctuations. 
Interestingly, a natural distinction appears between a normal regime, where 
price fluctuations are moderate, from a 'crash' regime, where panic effects 
are self-reinforcing and leads to a price collapse. 



2 A Short Survey of Empirical Data 
2.1 Linear correlations 

We shall denote in the following the value of the stock (or any other asset) 
at time t as x(t), and the variation of the stock on a given time interval 
r as 5 T x(t) = x(t + r) — x(t). The time delay r can be as small as a few 
seconds in actively traded markets. However, on these short time scales, 
the fluctuations cannot be considered to be independent. For example, the 
second order correlation function defined as Q: 

is significantly non-zero up to At = 15 minutes, as shown in Fig. 1. 

Therefore, it is certainly inadequate to think of the price process as a 
'martingale' for small time increments. This however does not necessarily 
mean that there are arbitrage opportunities: one can easily see that even 
very small transaction costs prevent the use of these short time correlations, 



^dn principle, the average value of 5 T x(t) should be removed, but this leads to in- 
significant corrections on the time scales considered. For the same reason, we neglect the 
difference between the variation of the price and the more often studied variation of the 
log of the price. 
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Figure 1: Correlation function C(At) of the one minute price increments of the 
S&P 500, for the period 1983-1996, which decays to zero on the scale of r* ~ 15 
minutes. Note that r* has drifted to smaller values with time. 
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which would imply a very high (and thus very costly) trading frequency M. 
For time delays At longer than a certain r*, say one hour, the autocorrelation 
of price increments is very nearly zero; correspondingly, the variance of price 
increments grow linearly with At for At > t*. 



2.2 Distribution of elementary increments 

The simplest idea is thus that the increments f\ 5(t) are independent identi- 
cally distributed (hd) random variables. In this case, the knowledge of the 
distribution density P* of 5 would suffice to reconstruct the distribution of 
increments on any time delay At larger than t*, through a simple convolu- 
tion. The shape of P* is strongly non gaussian: an estimate of its kurtosis 
K T * on a historical basis leads to numbers on the order of 20 (again on liquid 
markets). As shown in Fig. 2 on the example of the British Pound/U.S. $ 
time series, the tails of P* decay as an exponential: 

P(5) oc exp(-A±|5|) 5^±oo (2) 

(although several authors have reported a somewhat slower decay, as a power 
law l^l -1 ^ with a rather large exponent v ~ 3-4 |7[], or possibly a 'stretched 
exponential' M). Such slowly decaying tails survive upon convolution, and 



are thus particularly relevant for extreme risks forecasts JTO . 

A reasonable fit of P* on most markets can be achieved using a symmet- 
rical 'truncated' Levy distribution ||, ||, defined in Fourier space as: 

\ogP*(z) = A \\» - (A 2 + z 2 r /2 /Wl ; /(*) = cos (//arctan M ] 

COS7T/i/2 L J \ A J 

(3) 

where A the scale factor, \i is the 'tail exponent', which is found to be close to 
1.5 for all markets, and A describes the exponential fall of the far-tails. Note 
that when \x = 2, one recovers the gaussian distribution, while for A — > 0, 
one finds the stable Levy distributions proposed by Mandelbrot with tails 
which decay as |5| _1_At . 

2 In the following, we shall simplify the notation and use 6(t) = S r *x(t) for the incre- 
ments on the time scale t*. 



5 



Distribution of returns (BP/$) 

Exponential + Power-law fit 



10000 



1000 



100 



o 



10 




1 

0.000 



0.010 



0.020 
Amplitude 



0.030 



0.040 



Figure 2: Distribution of the British Pound/US $ negative increments for r* = 
30 mimutes, in the period 1991-1995. The 'near-tail' of the distribution can be 
fitted by a power-law, which clearly overestimates the 'far-tail'. The latter is well 
represented by an exponential fall-off, Eq. d2T): note the linear- log scale 
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2.3 Anomalous decay of the kurtosis and volatility per- 
sistence 

The N th autoconvolution of P*, where N = At/r* is simply obtained by 
scaling the parameter A to NA. With no further ajustable parameters, 
this leads to a fair representation of the historical distribution of S^t^ [§]; 
|TT|| . However, some systematic differences (in particular in the tails of the 
distribution) show up, which reveal the inadequacy of such a simple iid 
hypothesis. For example, it is easy to show that under such an hypothesis, 
the kurtosis of the increments on scale At = Nt* should be a factor iV 
smaller than the kurtosis on scale r*. Empirically, however, one finds a much 
slower decay, as kn = k>i/N x , with A in the range 0.3 — 0.6 (see Fig. 3). This 
can be related to another observation, which was made many times in the 
literature || [12]: higher order correlation functions of the price increments, 



such as: 

{5\t + mr*)5\t)) 

= ^ 1 ( 4 ) 

decay only very slowly with time. A simple fit of this decay is again of the 
power law type |5], [13], 

C 2 (m) ~ (5) 

with a exponent A' in the same range as A. A simple way to rationalize these 
findings is to assume that while the sign of the price increment is completely 
decorrelated as soon as At > r*, the amplitude of the increment (which is 
a measure of the market activity) is correlated in time. Bursts of market 
activity, related to external news or crisis, can easily persist for several days, 
sometimes months - this explains why C^m) decays rather slowly. We shall 
thus assume that the increment 5 can be written as: 

6(U) = 7(*t) x e{U) U = it* (6) 

where the scale 7 > measures the amplitude of the increment, which can 
be thought of as the local volatility of the market. The random variable e is 
short range correlated (over time r*) of mean zero and variance unity (and 
independent from 7). Note that e is not necessarily Gaussian, as assumed in 
ARCH-like models ||12||. It is then rather easy to show that the kurtosis of the 
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increments on scale At = Nt* can be expressed as [IR 



K N 



N 



N 



m 

n + (3 + Ko )C 2 (0) + 6 £ (1 - -)C 2 (m) 

m=l iV 



(7) 



where k is the kurtosis of the random variable e. Note that even if k = 0, 
the kurtosis of P* is non zero due to the randomly fluctuating scale parameter 
7. If one furthermore assumes that C2 decays as a power-law (Eq. |5|), then 
from Eq. (|7|), k n decays for large N as N~ x with A = A'. 



2.4 Apparent (multi)-scaling behaviour 

Hence, the fact that the scale 7 of the random increments has long range 
temporal correlations induces an anomalously slow convergence of the sum 
x(Nt*) — x(0) = J2iLi$(U) towards the gaussian distribution. One should 
note that this can induce apparent scaling behaviour on restricted time inter- 
vals. For example, a numerical simulation of the random scale model (||) with 
a slowly decaying C^m) can be analyzed in scaling terms [IB , i.e. fitting 
the moments of x(Nt*) — x(0) as power laws: 

(\x(Nt*) -x(0)\ q ) oc#« (8) 

This actually works quite well, and leads to a family of exponents Q which 
deviates from the theoretical straight line ( q = q/2 which holds in the limit 
iV — > 00, provided all the moments of 7 exist. For q = 2, the relation C2 = 1 
holds because the linear correlation function C(At) is zero for At > r*. For 
q = 4 one finds, using the definition of the kurtosis: 

([x(Nt*) - x(0)] 4 ) = (3 + k n )([x(Nt*) - x(0)] 2 ) 2 oc 3iV 2 + k x N 2 - x (9) 

which is thus the sum of two power laws. This can be however fitted with a 
unique 'effective' exponent £4, which can be substantially below the asymp- 
totic value = 2 since N~ x is not very small in practice. This argument 
holds true for higher moments for which ( q < q/2; the difference between 
( q and q/2 actually grows with q (for a fixed value of N). We believe that 
this might be the reason for the 'multifractal' scaling recently reported in the 
literature [117, 
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2.5 Time reversal symmetry 



The last empirical fact which we would like to comment on is the question 
of time reversal symmetry. Is it possible to detect the direction of time 
in a financial time series ? The answer to this question would be no for 
a simple Brownian motion, for example, but also for a much wider class 
of processes, such as the multifractal time construction proposed in [|Hf|. 
Correlation functions sensitive to the arrow of time have been proposed by 
Pomeau |19|]. One example is: 



C T (At) = (x{t) [x(t + At) - x{t + 2At)} x{t + 3 At)) (10) 

which is non zero if the time triplets t,t + At, t + 3At and t,t + 2At, t + 3 At 
cannot be distinguished statistically. For the price series itself, the above 
correlation function is zero within error bars. But when one studies the 
'volatility' process 7, then this skew correlation function is distinctively non 
zero, and reaches a maximum (in the case of the S&P 500) for At ~ 1 
month. This shows that the volatility time series is not invariant under time 



reversal. A similar conclusion is also reported in [20 , 14], where the authors 
observe that a high 'coarse-grained' volatility in the past increases in a causal 
way today's 'fine-grained' volatility. This is not unreasonable, as one feels 
intuitively that an anomalously large change of the close to close price over - 
say - the previous week triggers more intraday activity the following week. In 
any case, we feel that this absence of time reversal symmetry is an important 
(albeit less frequently discussed) stylized fact of financial time series. 



3 Implications for Option Pricing 

3.1 General framework 

We now turn to the problem of option pricing and hedging when the statis- 
tics for price increments have the non-Gaussian properties discussed above. 
The distinctive feature of the continuous time random walk model usually 
considered in the theory of option pricing is the possibility of perfect hedging 



[pT[| , that is, a complete elimination of the risk associated to option trading 



This property however no longer holds for more realistic models [22 . 
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Let us write down the global wealth balance AH^Iq associated with the 
writing of a 'call' option of maturity T and exercice price x s @: 

AW\% = C{x ,x s ,T)exp{rT) - max(x(T) - x s ,0) (11) 
+ tfatJexviriT-tiWi-rxitjT], (12) 

i 

where C(xo,x s ,T) is the price of the call, xo = x(t = 0) and <f>(x,t) the 
trading strategy, i.e. the number of stocks per option in the portfolio of 
the option writer. Finally, r the (constant) interest rate. The second term 
defines the option contract: the profit of the buyer of the option is equal to 
x s — x(T) if x(T) > x s (i.e. if the option is exercised) and zero otherwise: the 
option is an insurance contract which guarantees to its owner a maximum 
price for acquiring a certain stock at time T. Conversely, a 'put' option would 
guarantee a certain minimum price for the stock held by the owner of the 
option. 

A natural procedure to fix the price of the option C(xq, x s ,T) and the 
optimal strategy <p*(x,t) was proposed independently in |23|, |22[ and further 
discussed in || p5| . It consists in imposing a fair game condition, i.e.: 

(AW\M) = (13) 

and a risk minimisation condition: 

5(AW\^} 



8(j)(x, t) 



(14) 



Here, we assume that the variance of the wealth variation is a relevant mea- 
sure of the risk. However, other measures are possible, in particular the 
'Value-at-Risk', which is directly related to the weight contained in the neg- 
ative tails of the distribution of AVT|J. 



The notation (...) in Eqs. ( JTS] , ^ ) means that one averages over the prob- 
ability of the different trajectories. The explicit solution of Eqs. (0,0) for 
a general uncorrelated process (i.e. (SiSj) = for i ^ j) is relatively easy to 
write if the average bias (5) and the interest rate r are negligible^, which is 
the case for short maturities T. In this case, one finds 



C(x ,x s ,T) = / dx' (x' -x s )P(x',T\x o ,0) (15) 

J X„ 



3 For the general case, see (q, Eil 
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POO ( rp' T \ 

<P*(x,t) = dx'{8) {x , t) ^ {x , }T) K °> P{x',T\x,t) (16) 

Jx 3 0~ yX, t J 

where o~ 2 (x, t) = {b~ 2 )\ x> t is the 'local volatility' - which may depend on x,t 

- and (S)( x ,t) >(x',t) is the mean instantaneous increment conditioned to the 

initial condition (x,t) and a final condition (x',T). The minimal residual 
risk, defined as TV = (AW\q [4>*] 2 ) is in general strictly positive (and in 
practice rather large), except for Gaussian fluctuations in the continuous 
limit, where the residual risk is strictly zero! In this limit, the above equations 
T3| , |H| ) actually exactly lead to the celebrated Black-Scholes option pricing 



formula. In particular, one can check that 0* is related to C through: 
dC(x ,x s ,T)/dx . 



3.2 Cumulant expansion and volatility smile 

In the case where the market fluctuations are moderately non-Gaussian, one 
might expect that a cumulant expansion around the Black-Scholes formula 
leads to interesting results. This cumulant expansion has been worked out 
in general in 0. If one only retains the leading order correction which is 
(for symmetric fluctuations) proportional to the kurtosis, one finds that the 
price of options C(xo,x s ,T) can be written as a Black-Scholes formula, but 
with a modified value of the volatility a, which becomes price and maturity 
dependent |15| : 




(17) 



The volatility o~i mp . is called the implied volatility by the market operators, 
who use the standard Black-Scholes formula to price options, but with a 
value of the volatility which they estimate intuitively, and which turns out 
to depend on the exercice price in a roughly parabolic manner, as indeed 
suggested by Eq. (|17|). 

This is the so-called 'volatility smile'. Eq. flT7|) furthermore shows that 
the curvature of the smile is directly related to the kurtosis of the underlying 
statistical process on the scale of the maturity T = Nt*. We have tested this 
prediction by directly comparing the 'implied kurtosis', obtained by extract- 
ing from real option prices (on the BUND market) the volatility a (which 
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Figure 3: Plot (in log-log coordinates) of the average implied kurtosis Ki mp (de- 
termined by fitting the implied volatility for a fixed maturity by a parabola) and 
of the empirical kurtosis kn (determined directly from the historical movements 
of the BUND contract), as a function of the reduced time scale iV = T/r, r = 30 
minutes. All transactions of options on the BUND future from 1993 to 1995 were 
analyzed along with 5 minute tick data of the bund future for the same period. We 
show for comparison a fit with kn — N~ - 6 (dark line). A fit with an exponentially 
decaying C2(m) is however also acceptable (dotted line). 
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turns out to be highly correlated with a short time filter of the historical 
volatility), and the curvature of the implied volatility smile, to the historical 
value of the kurtosis k n discussed above. The result is plotted in Fig. 3, 
with no further ajustable parameter. The remarkable agreement between 
the implied and historical kurtosis, and the fact that they evolve similarly 
with maturity, shows that the market as a whole is able to correct (by trial 
and errors) the inadequacies of the Black-Scholes formula, and to encode in 
a satisfactory way both the fact that the distribution has a positive kurto- 
sis, and that this kurtosis decays in an anomalous fashion due to volatility 
persistence effects. However, the real risks associated with option trading 
are, at present, not satisfactorily estimated. In particular, most risk control 
softwares dealing with option books are based on a Gaussian description of 
the fluctuations. 



4 A Langevin Approach to Price Fluctuations: 
Feedback Effects 

In the author's opinion, there is as yet no convincing 'microscopic' model 
which explains the distinctive statistical features of price fluctuations which 
were summarized in section [|, although many proposals have been put for- 
ward []25|, |27|] . In the spirit of statistical physics, it might be possible to 



describe market dynamics at a level of description which is intermediate 
between the macroeconomic level which is that of the market equilibrium 
models [f26j and the individual agent level which is that of the market mi- 
crostructure theory. Our basic idea is that although the modelling of each 
individual participants ('agents') is impossible in quantitative terms, the col- 
lective behavior of the market and its impact on the price in particular can 
be represented in statistical terms by a few number of terms in a (stochastic) 
dynamical equation. Our approach is in the spirit of many phenomenolog- 



ical, 'Landau-like' approaches to physical phenomena We have thus 

proposed a description of the dynamics of speculative markets with a simple 
Langevin equation [p9fl . This equation encapsulates what we believe to be 



some essential ingredients; in particular, the feedback of the price fluctua- 
tions on the behaviour of the market participants. The construction of a 
convincing model for large fluctuations is important for risk control, since 
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empirical data may be unreliable in the extreme tails of the distributions. 



4.1 A Phenomenological Langevin Equation 

At any given instant of time, there is a certain number of 'buyers' which we 
call 4>+(t) (the demand) and 'sellers', 4>-(t) (the supply). The first dynamical 
equation describes the effect of an offset between supply and demand, which 
tends to push the price up (if + > <p_) or down in the other case. In general, 
one can writef]: 

— = T(A<f>) A</> := 0+ - 0_ (18) 

where T is an increasing function, such that !F(0) = 0. In the following, we 
will frequently assume that T is linear (or else that A0 is small enough to 
be satisfied with the first term in the Taylor expansion of J 7 ), and write 

where A is a measure of market depth i.e. the excess demand required to 
move the price by one unit. When A is high, the market can 'absorb' sup- 
ply/demand offsets by very small price changes. Now, we try to construct 
a dynamical equation for the supply and demand separately. Consider for 
example the number of buyers <f> + . Between t and t + dt, a certain fraction of 
those get their deal and disappear (at least temporarily). This deal is usually 
ensured by market makers, which act as intermediaries between buyers and 
sellers. The role of market markers is to absorb the demand (and supply) 
even if these do not match perfectly. The effect of market makers (mm) can 
thus be modelled as: 

-r±0± (20) 



dt 

where T± are rates (inverse time scales). We furthermore assume that market 
makers act symmetrically, i.e, that T + = T_ = 7. On liquid markets, the 
time scale I/7 before which a deal is reached is short; typically a few minutes 
(see also below for another interpretation of I/7). 

There are several other effects which must be modeled to account for 
the time evolution of supply and demand. One is the spontaneous (sp) 

4 For simplicity, we use in the following a continuous time formalism, although discrete 
time evolution equations would be more adapted. 
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appearance of new buyers (or sellers), under the influence of new information, 
individual need for cash, or particular investment strategies. This can be 
modelled as a white noise term (not necessarily Gaussian): 



dt 



m±{t)+r]±{t) 



(21) 



where r)± have zero mean, and a short correlation time r c . m± is the average 
increase of demand (or supply), which might also depends on time through 
the time dependent anticipated return R(t) and the anticipated risk £(i). It 
is quite clear that both these quantities are constantly reestimated by the 
market participants, with a strong influence of the recent past. For example, 
'trend followers' extrapolate a local trend into the future. On the other hand, 
'fundamental analysts' estimate what they believe to be the 'true' price of 
the stock; if the observed price is above this 'true' price, the anticipated 
trend is reduced, and vice-versa. In mathematical terms, these effects can be 
represented as: 



R(t) = R + a f dt'K R (t - t 

J —oo 



, dx 



dt 



- - K{X - X ) 



(22) 



where K R is a certain kernel (of integral one) defining how the past average 
trend is measured by the agents, and k is a mean-reversion force, towards 
the average (over the fundamental analysts) 'true price' xq f\ 

Similarly, the anticipated risk has a short time scale contribution. It is 
well known that an increase of volatility is badly felt by the agents, who 
immediately increase their estimate of risk. Hence, we write: 



£(t) =S + /3 /* dt'Kx(t-t') 

J —oo 



dx 
dP 



Correspondingly, expanding m±(R, E) to lowest order, one has: 



/t djr ft 

dt'Kf>(t — t') — +j3± / dt K^(t—t) 
-oo dt J —oo 



dx 
dt 



(23) 



-K±(x-Xq) 

(24) 



5 Note that x is actually itself time dependent, although its evolution in general takes 
place over rather long time scales (years). 
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where the signs of the different coefficients are set by the observation that 
m + is an increasing function of return R and a decreasing function of risk 
E, and vice versa for m_. Eq. ( 23 ) contains the leading order terms which 
arise if one assumes that the agents try to reach a tradeoff between risk 
and return: the demand for an asset decreases if is recent evolution shows 
high volatility and increases if it shows an upward trend. This is the case 
for example if the investors follow a mean- variance optimisation scheme with 
adaptive estimates of risk and return ||, or the Black-Scholes option hedging 
strategy Actually, the 1987 crash is often attributed (at least in part) to 
the automatic use of the Black-Scholes hedging strategy, which automatically 
generates sell orders when the value of the stock goes down. 

We are now in position to write an equation for the supply/demand offset 
A0 by summing all these different contributions: 



dA(j) 
dt 



/ dx 
dt'K R (t-t') — 
-oo CLt 

b f dt'K s {t - t') 

J —oo 



dx 



k(x - X ) + T](t) 



with a,b,k > 0. Note in particular that b > reflects the fact that agents 
are risk averse, and that an increase of the local volatility always leads to 
negative contribution to A0. This feature will be crucial in the following. 
For definiteness, we will consider 77 to be gaussian and normalize it as: 



{r)(t) V (t')) = 2X 2 DS(t - t') 



(25) 



where D measure the susceptibility of the market to the random external 
shocks, typically the arrival of information. As discussed in section [2|, D 
should also depend on the recent history, reflecting the fact that an increase 
in volatility induces a stronger reactivity of the market to external news. In 
the same spirit as above, one could thus write: 



D = D + Dt f dt'K D {t - t') 

J —oo 



dx 
dP 



(26) 



For simplicity, we will neglect the influence of D\ in the following. But such 
a term might be responsible for the time reversal symmetry breaking effect 
reported in section |2|, as well as for the fact that volatility appears to be time 
dependent. 
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4.2 Simple consequences. Liquid vs. Illiquid markets 

Let us consider the linear case 'risk neutral' case where 6 = 0. We will 
assume for simplicity that K^{t) = Texp— (rt), and first consider the local 
limit where T is much larger that 7 (short memory time). In this case, the 
equation for x becomes that of an harmonic oscillator: 

d 2 x ( a\ dx k , _ , 1 



*> + I 7 " A ) Tt + X {X ~ = A" W (27) 
where m has been absorbed into a redefinition of Xq := XQ+m Q /k. For liquid 
markets, where A and 7 are large enough, the 'friction' term 7 := 7 — a/A 
is positive. In this case the market is stable, and the price oscillates around 
an equilibrium value Xq, which is higher than the average fundamental price 
if the spontaneous demand is larger than the spontaneous supply (i.e. m 
is positive), as expected when the overall economy grows. The situation 
is rather different for illiquid markets, or when trend following effects are 
large, since 7 can be negative. In this case, the market is unstable, with an 
exponential rise or decay of the stock value, corresponding to a speculative 
bubble. However, in this case, dx/dt grows with time and it soon becomes 
untenable to neglect the higher order terms, in particular the risk aversion 
term proportional to b, which is responsible for a sudden collapse p9 . 



4.3 Risk aversion induced crashes 

In the case of liquid markets on short time scales, it is reasonnable to set 
k = PPI. Setting u = dx/dt and still focusing on the limit where the 



memory time T 1 is very small, one finds the following non linear Langevin 
equation: 

du m _ b 9 1 , . dV 1 . . 

_ = _£_ 7u __„ 2 + _ I)(t ),__ + _, (t ) (28) 

This equation represents the evolution of the position u of a viscous ficti- 
tious particle in a 'potential' V(u). In order to keep the mathematical form 
simple, we set the average trend mo/ A to zero (no net average offset between 
spontaneous demand and spontaneous supply); this does not qualitatively 
change the following picture, unless m is negative and large. The potential 
V(u) can then be written as: 

V{u) = \u 2 + ^-u 3 (29) 
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which has a local minimum for u = 0, and a local maximum for u* = — A7/6, 
beyond which the potential plumets to —00. The 'barrier height' V* sep- 
arating the stable region around u = from the unstable region is given 



The nature of the motion of u in such a potential is the following: starting 
at u — 0, the particle has a random harmonic-like motion in the vicinity of 
u = until an 'activated' event (i.e. driven by the noise term) brings the 
particle near u*. Once this barrier is crossed, the fictitious particle reaches 
— 00 in finite time. In financial terms, the regime where u oscillates around 
u = and where b can be neglected, is the 'normal' random walk regime 
discussed in the previous paragraph. (Note that the random walk is biased 
when m 7^ 0). This normal regime can however be interrupted by 'crashes', 
where the time derivative of the price becomes very large and negative, due 
to the risk aversion term b which enhances the drop in the price. The point 
is that these two regimes can be clearly separated since the average time t* 
needed for such crashes to occur can be exponentially long, since it is given 
by the classical Arrhenius- Kramers formula [[301 , |3~T] : 



where D is the variance of the noise 77 and t\ = I/7. Taking t* = 10 years, 
a = 1% per day, and t\ = 10 minutes, one finds that the characteristic value 
u* beyond which the market 'panics' and where a crash situation appears is 
of the order of —1% in ten minutes, which not unreasonnable. Note that in 
this line of thought, a crash occurs because of an improbable succession of 
unfavorable events, and not due to a single large event in particular. Fur- 
thermore, there are no 'precursors' - characteristic patterns observed before 
the crash: before u has reached u*, it is impossible to decide whether it will 
do so or whether it will quietly come back in the 'normal' region u ~ 0. 
Note finally that an increase in the liquidity factor 7 reduces the probability 
of crashes. This is related to the stabilizing role of market makers, which 
appears very clearly. 

The above calculations can be extended to the case where the 'integration 
time' of the market (appearing in the kernels Kr and K-%) is not very short, 
with very similar conclusions [p5j] . 



by: 



V* = V(u*) - V(0) 




(30) 




(31) 



18 



5 Conclusion 



Research on financial markets can focus on rather different aspects. We have 
discussed here three interrelated themes: the statistical nature of the market 
fluctuations, its use for option pricing and the need for a microsocopic 'ex- 
plicative' model which accounts for empirical observations. We have empha- 
sized the fact that a good model of these fluctuations was crucial to estimate 
and control large risks. This requires in particular an adapted theory for op- 
tion pricing, which goes beyond the traditionnal Black-Scholes dogma. We 
have shown that subtle statistical effects, such as the persistence of volatility 
fluctuations, is rather well reflected in option prices. This shows that the 
market as a whole behaves as an adaptive system, able to correct (through 
trial and errors) the theoretical inadequacies of the Black-Scholes formula. 

The theoretical understanding of the tails of the distributions is further- 
more of fundamental importance, because by definition the empirical obser- 
vation of rare events leads does not allow one to determine the extreme tails 
of the distribution with great accuracy. In this respect, an interesting as- 
pect of the Langevin equation discussed above is that crashes events appear 
as fundamentally distinct from 'normal' events, and their probability of oc- 
curence is thus not expected to lie on the extrapolation of the (non Gaussian) 
distribution constructed from these 'normal' events. This is actually what 
is observed: although the empirical tail of the distribution of the S&P 500's 
daily increments is pretty well fitted by an exponential down to probabilities 
of 10~ 3 , events which have a probability of 10~ 4 (i.e. crashes which occcur 
once every 40 years) have an amplitude which is much larger than expected 

There are many other aspects of financial markets which are worth in- 
vestigating using the tools and ideas of statistical physics, which we have 
not discussed here at all: we refer the reader to the recent literature, which 



can be found on the preprint server "cond-mat" [33]. One of the most excit 



ing subject is the modelling of the interest curve, which might have a lot in 



common with the dynamics of an elastic string driven by noise [33|. 
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